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ABSTRACT  fConflnif*  ofi  r*v»r««  oido  it  n«c««««ry  «nd  Idoniity  by  block  numbor) 

Hiis  report  describes  four  receivers  for  the  detection  of  a sinusoid  of  known  frei)uency  Inii  unknown  phase  in 
(iaiissiun  white  noise  of  known  spectmm  level.  These  are;  the  likelihood  ratio  detector,  a puKessoi  emplox  ing 
mcoherenlly  averaged  Discrete  Fourier  Transforms  (DF  I s)  of  the  received  signal,  and  two  systems  utili/ing  the 
fre»|nency  response  of  adaptive  1 meai  Ihedictive  Filteis  (1 1’Fs).  11ie  I I'F  systems  ate  desciibed  in  detail,  and  an 
expression  for  optimal  adaptive  (liter  parameter  .selection  is  derived.  Power  curves  for  the  four  receivers  are 
derived  and  compared.  Results  indicate  that  the  adaptive  methods  are  similai  in  peiformance  to  the  mcoheient  DF  l 
proces-sor  with  one  of  the  adaptive  techniques  offering  a I dll  reduction  in  input  signal  lonoise  ratio  (.SNR)  foi 

_ ‘ M lIL 
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I.  INTRODUCTION 


The  use  of  an  adaptive  linear  predictive  filter  to  detect  the  presence  of  a sinusoid  in 
white  noise  has  been  discussed  by  Widrow  t7  ul  (Ref.  1 ) in  a comparison  of  two  detection 
methods,  I he  first  uses  as  its  test  statistic  the  frequency  response  magnitude  of  an  adaptive 
linear  predictive  filter  (LIM’),  the  second  employs  incoherently  averaged  discrete  Fourier 
transforms  (DF'Ts)  of  the  input  data.  Widrow's  comparison  is  based  on  output  signal-to- 
noise  figures-of-merit  (FOM)  for  the  two  processors  and  indicates  that  the  LI’F  technique 
may  offer  significant  improvement  over  the  DFT  when  incoherent  averaging  is  required  to 
process  the  available  data. 

I ufts  (Ref.  2)  has  challenged  Widrow’s  analysis  on  the  grounds  that  the  DIT  should 
have  been  allowed  to  process  the  data  coherently.  In  their  reply.  Widrow  ci  al  (Ref.  .^)  have 

agreed  that  a detector  utilizing  a coherent  DFT  of  the  data  would  indeed  perform  better  ' 

than  the  Li’F  detector  and  suggested  that  a different  FOM  is  appropriate  for  describing  I.PI'  | 

performance.  i 

Despite  this  interchange,  the  fundamental  question  of  detector  performance  remains 
unresolved.  .Analysis  published  to  date  (Ref.  I ) is  based  on  FOM  arguments,  and  these  can 
be  misleading  when  used  in  the  comparison  of  processors  with  different  statistical  properties. 

In  order  to  provide  a definitive  answer,  the  problem  should  be  addressed  from  the  standpoint 
of  decision  theory  (Refs.  4,.‘i ).  This  paper  represents  a further  step  toward  defining  the  per- 
formance of  detectors  utilizing  LPFs  by  first  c.xtending  Widrow’s  analysis  to  account  for 
adaptation  of  (he  LPF,  then  using  statistical  decision  theory  to  detennine  pcrl'ormancc. 


II.  PROBLEM  DESCRIPTION 

riie  detection  of  a sinusoid  in  Ciaussian  white  noise  is  to  be  considered.  I'he  fretiuency 
of  the  sinusoid,  fj,,  total  noise  power,  , and  noise  power  per  hertz,  o? , are  known.  Phase 
and  amplitude  of  the  sinusoid  are  unknown,  but  are  constant  over  the  observation  interval. 
Fach  detector  is  to  observe  the  data  over  a fixed  interval,  then  decide  which  of  the  two 
possible  hypotheses  is  true: 

llo  noise  alone  present 

H,  sine  wave  and  noLse  present. 

The  a priori  probabilities  of  these  hypotheses  are  unknown.  Hence,  a Neyman- 
Pearson  test  (Ref.  4)  will  be  used  to  decide  whether  Hq  or  II,  is  true.  A scalar  test  statistic 
generated  by  each  detector  will  be  compared  to  a fixed  threshold  for  that  detector.  The 
threshold  is  set  to  achieve  a desired  probability  of  false  alarm  (i  t .,  the  probability  ol  choos- 


l our  iloti'ctor  siructurcs  tlcscrihed  in  the  lollowiiij!  section  ;ire  to  be  eonsuleieii 

1 1 ) The  optimal  lieteelor  (I  ig.  1 1 

(2)  A\er.igcd  spectral  estimates  ol  the  data  (l-ig.  2). 

t3)  Fourier  translorm  ol  the  LI’F  weight  vector  (Fig.  31. 

(4)  Fourier  transform  magnitude  sipiared  of  the  1 1*1'  weight  vector  (Fig.  4). 

llie  predictive  filters  of  Figs.  3 and  4 are  described  in  the  Appendix. 
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Figure  1 . Coherent  DFT  detector  for  a sinusoid  of 
frequency  f^  in  bandliiniled  Caussian  white  noise  of 
known  s(>ectruin  level. 
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sinusoid  of  tiequenev  f^  in  handlinuted  Ciaussian  noise  of 
known  spectiuin  level 


I 


1 


SAMPLE  At 

t a 


I 1,2,.  .N 


INPUT 


I DELAY  j 

, r- 

I samples] 

FILTER 

WEIGHT  VECTOR 


ADAPTIVE 

FILTER 


WEIGHT 

UPDATE 

ALGORITHM 


[dpt] 


DFT  AT 

FREQUENCY  f, 

; PHASE  ROTATION 
i exp  [-)2:r1,A<S| 

REAL  PART 

: ri-:: 

-y  ■ Tj  COMPARATOR 
I [^THRESHOLD  T 

Ho 


h" 


Hi 


SAMPLE  AT 


lA 

1.  2. 


..  N 


INPUT 


DELAY 

T 

SAMPLES 


ADAPTIVE 

FILTER 


FILTER 

WEIGHT  VECTOR 


WEIGHT 

UPDATE 

ALGORITHM 


DFT  AT 

FREQUENCY  l,| 

I MAGNITUDE  SQUARe] 

y tJ  comparator  It  t 

j[  THRESHOLD  - TJ^ 
Hn  H, 


Figure  .1.  Coherent  linear  predictor  (CLP)  detector  Figure  4.  Incoherent  linear  predictor  (ILP)  detector 

for  a sinusoid  of  frequency  f^  in  bandliiniled  for  a sinusoid  of  frequency  f^  in  bandliniited 

Gaussian  noise  of  known  sirectruni  level.  Gaussian  noise  of  known  spectrum  level. 


111.  ANALYSIS 


OPTIMAL  DETECTOR  - COHERENT  DFT 

Titc  likelihood  ratio  detector  for  a signal  of  unknown  amplitude  and  phase  is 
discussed  by  Helstrom  (Ref.  .^).  Tliis  optimal  detector,  which  is  the  coherent  DFT  mentioned 
by  Tufts  ( Ref.  2).  can  be  implemented  by  splitting  the  input  data  into  two  channels  as  shown 
in  F'ig.  1 . The  first  channel  is  multiplied  by  cos(27rfj,t)  and  integrated  over  the  observation 
interval.  The  second  channel  is  multiplied  by  sin(2rrfj,t)  and  integrated.  At  the  end  of  the 
observation  period  a test  statistic,  denoted  by  7('oH  P'S-  R>rmed  by  squaring  and 

summing  the  two  integrals.  ITic  Neyman-Pearson  test  is  applied  to  7('qu  in  order  to  decide 
signal  presence  or  absence. 

One  practical  method  of  implementing  this  detector  uses  the  DF  T.  lire  observed  data 
are  sampled  at  Nyquist  rate  and  the  msulting  data  .sequence  is  transformed  with  the  known 
sinusoid  frequency,  fj..  at  one  of  the  DFT  frequencies.  Real  and  imaginary  parts  of  the  DFT 
at  this  frequency  are  squared  and  summed  to  form  the  test  statistic.  Note  that  the  data  must 
be  transfornted  coherently.  Tliat  is.  either  a single  DFT  must  process  the  entire  data  sequence 
in  one  transformation,  or  short  DFTscan  be  summed  in  phase  before  the  test  statistic  is 
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I OHERtSr  FOURIER  TRANSFORM  OF  THE  LPF  WEIGHT  VEC  TOR  (C  LP> 


IIk'  I IM  used  in  this  work  is  discussed  in  the  Appendix  llie  weight  \ecior  is 
modeled  as  a Miuisisivl  oT  lmie-\aryin^  amplitude  summed  with  a Gaussian  “misad.iustment" 
noise  1 4|uations  (22).  (2.^).  (25).  and  ( 2(')  describe  the  TIM'  weiuhl  vector.  Wti.k).  where 
I denotes  tlie  individual  weight  number  (0  ^ i < N-l ) and  k denotes  the  number  ol  adapta- 
tion iterations  tk  = 1 ,2.  • • • ).  The  C'l  P iletector  (Tig.  .))  uses  as  its  test  statistic  the 
measured  amplitude  ot  the  weight  vector  cosine  component  This  amplitiule  is  estimaied 
by  the  sum 


N-l 

LI’  “ S Wti.k)cosl2»rl'j;ti  ♦ 6)  A). 
1-0 


t(i) 


where  is  the  sinusoid  rrequency  to  be  detected.  5 is  tlie  l.l’l'  prediction  delay  tl  ig  .G.  and 
A IS  the  sampling  period  of  the  digital  system.  The  test  statistic  is  computed,  as  illustrated 
in  Tig.  .V  by  first  forming  the  OFT  of  the  weight  vector  ti  e.,  evaluating  the  l.l’l'  frequency 
response  at  frequency  fj,).  then  rotating  the  phase  of  this  complex  quantity  as  indicated. 

Tlie  real  part  of  the  result  is  used  as  the  test  statistic.  Since  the  weight  vector  misadjiistment 
noise  is  assumed  to  be  Gaussian,  the  test  statistic  will  also  be  Gaus.sian.  llie  variance  of  the 
misadjustment.  and  hence  the  variance  of  the  test  statistic,  is  constant  under  Hj,  and  H, . 

In  the  present  analysis  it  is  assumed  that  all  weights  are  initially  set  to  zero  ti  e.,  NV’ti.O)  = 

0.  0 < i < N-l ) and  that  a single  OFT  is  calculated  at  the  end  ot  the  obseixation  period 
When  nonnali/ed  to  unit  variance,  the  probability  distributions  of  the  test  statistic  under 


H, 

are  (see  .Appendix) 
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P-y^.ll,  tG  = Nfla.l). 

tS) 

where  Nt  • x.y)  denotes  the  normal  density  with  mean  x and  variance  y. 

N N,- 

1 - 11  - :no-(  \ SNR)1  * 


f N 2|  'j 
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and  M the  adaptive  feedback  constant  of  the  I I’F'. 


N 

— SNR 
N 

1 +-  SNR 
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llie  probability  of  sinusoid  defection  can  Ih'  maximized  bv  si'paiating  the  test 
statistic  distributions  under  Hq  and  ll|  as  far  as  possible,  i c . choosing  the  controllable  1 I I’ 
parameters  \ and  n so  as  to  make  o as  large  as  possible  .Assuming  N to  Iv  tixed  bv  haulwaie 
limitations,  this  is  accomplished  by  choosing 


1.25^4.^ 
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SNK  IN  a uoiumal  value  v'l  SNK  tv'  be  iUMoeU\i.  auJ  tlie  uumeiatv'i  ol  I q < 101  is  llu' 
iu>ii/eii>  u>i>t  v'l  Ihe  lianseeiulemal  ei|uatuM\  1 1 + J\W"'  1 \Vul\  llus  eboiee  ot  m- 

vlolev'lv'i  I'l'tlonuaiKc  i>  v'ptmu/cvl  when  SNR  S!^R  I'uiveN  pieNeiuevl  in  Seelion  l\  will 
show  that  Oeteelor  pertv'inianee  in  not  Ntuvnislv  inllueneevl  bv  sSiR.  aiul  Iheielv'iv  itN  value 
IN  not  entieal 

SubNtitutin^  K)  1 101  intv'  1\|n  t'>l  aiul  tSl  jsiveN  the  ileuMtv  iM  the  test  NtatiNtie  aN  a 
tiinetion  ol  N.  N | . SNR.  aiul  SNR  lluAe  parainetetN  vletermine  the  t I T vleteetoi 
|Vi  t'orinanee. 


IMOHl  RFNI  KOl'RlbR  IRANSI  ORMOl  1 III  I PI  Wl  U'.Ml  VI  l lOR  (II Pi 

Hus  ileteetor  d m 41  Ows  not  utilize  the  deteiniintstie  phase  ol  the  1.1*1  vvei.cht 
veetor  I'he  test  statistie  is  the  nucnitiule  siiuareO  ol  the  vvemhi  veelv'r  1>1'  T (i  e . the 
I PI'  livviueney  rv'sponse  inagnituvle  snuaiwll  at  tieiiuenev  1,^  llie  model  I'l  the  wemhl 
veetor  IS  develo(vd  in  the  .Appendix.  I'robabihty  densities  ol  the  test  statistie  under 
hypotheses  Hj,  and  II,  are  (with  normalization  to  unit  varianeel 


tin 


ti:i 


where  a is  jtiven  by  I q t*>V  .lust  as  lor  the  (.'1 1*  deteetv'i  . prv'bability  of  sinusoid  deteetion 
IS  maximized  by  elu'v'smv;  the  .ulaptive  feedb.iek  eonstant  aeev'rdinj:  to  I'q  (lOV 


IV.  PFRI  ORM  ANI  I f OMP ARISONS 

The  density  funetions  of  Seetion  111  will  now  be  used  to  develop  eomparisv'us  of 
deteetor  performanee  Results  are  based  on  the  Ne>  man-Peaixoii  hvpolhe.sis  test  and  aiv  m 
the  form  of  power  eurves  slu'vvinjs  deteetion  probabilities  as  a funv  tion  v'f  false  alarm  piob.i 
bihtv  and  SNR  of  the  reeeived  signal.  Ilie  plots  published  by  Mareum  tRef.  t'l  and  the 
nomogram  of  I'rkiwvitz  tRef.  '1  were  usevl  in  obtaining  these  eurves. 

l igure  5 ev'inpares  the  ev'hereni  Ol'  1,  ASP  , and  11 1*  ileteetvMs  for  the  ease  v'f 

10.2. 4tHi  input  data  Namples.  Hie  ASP  OP  P length  v'f  10*4  pv'ints  w.in  ehosen  to  evjual  the 

X 

I PI  weight  veetv'i  length  Hie  parameter  SNR  reviuired  m v'lder  to  elu'v'se  n toi  the  11  P 
w.iN  seleeled  evpial  to  1 S 10'-^  t-.Hl  dlO  l iguie  5 shows  that  the  11  P and  ASP  piv'eesNOiN 
ate  rv'Ughlv  evpial  in  performanee  l ive  11  P offers  a fraelu'ii  of  a vlM  .idvant.ige  at  the  higher 
SNR  values  shown,  but  the  ASl  peifornis  shghtlv  better  at  smaller  SNRs  Oils  result  diffeis 
ftv'in  the  ev'iiehisu'ii  v'l  NNulrv'w's  analvsiN  tRef  1 1 based  v'li  a 1 t>M  eonvp.iiisv'n  vvhieh 
indiv .liv'd  a large  diflereuv'e  between  the  kSl  and  11  P pu'eessv'is  .As  piedieled  bv  1 lifts 
(Ret  Jl.  the  ev'hereni  OP" P deteetor  in  markedlv  bettei  th.m  either  the  ASP  or  11  P delv'eiv'is 


COHERENT 


INPUT  SNR  (ilB) 

Tigurc  >.  I’onipuiison  oTcohoroiit  ITET,  11  P,  and  ASl 
detectoi  {■KriToiinance. 


f igurc  (1  complies  ll>c  t;ohcri’nt  Ol  T.  ASl'.  aiui  C'l  I’  lU'tcctors  tor  tlic  c.isc  of 
ID-. 400  input  Jafa  sjinplcs.  .M)  paraiDi'U'rs  arc  thi'  s;inu'  as  those  of  l-i};.  5.  I'iijurc  0 shovss 
ttial  lire  C l.P  offers  about  a l-ilB  mrpro\eineirt  oxer  tire  .AS!:  proeessor.  but  is  still  far  sliort 
of  optiirral  perforinaiiec. 

A 

I 1 sure  ' illustrates  the  effeet  of  varyiiis  SNR  on  (’l  l’  perfornianee.  lire  ehanse  in 
A ^ I 

IH'rfornranee  is  neglisible  as  SNR  is  reilueeil  from  10“'  to  10  . lire  SNR  renuireil  to  aehiexe 

\ A 

a si'x'U  I’d  in  tire  ranee  illustrated  inereases  onl\  about  dR  as  SNR  is  inereased  from  10“'  to 
10  “ Thus  the  ehoiee  ot  SNR  is  not  eritieal  to  I I’l  deteetor  pertormanee  in  this  example 

f igure  S shows  the  effeet  of  varying  the  weight  xeetor  length  of  the  C’l !’  deteetor  and 
PIT  length  of  the  ASl  deteetor.  The  relative  ivrformlTNee  of  the  two  deteetors  is  not 
sigmfieantlx  affeeted.  In  both  eases  the  ( I I’  deteetor  offers  about  a InlH  advantage  in  the 
SNR  required  to  aelueve  a given  I’p. 

f igure  eompares  .ASl  and  C’l !’  deteetor  pertormanee  for  two  values  of  N p,  the 
total  number  of  data  srimples  proeessed.  Relative  pertormanee  of  the  two  processors  is  not 
• strongly  affeeted.  At  tlie  smaller  value  of  N p.  the  ( I I'  deteetor  offers  shghtlv  more  than  a 
1-dB  advantage  in  the  SNR  required  to  aehieve  a given  •’iv 
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Figure  6.  Comparison  of  coherent  DFT,  ASF. 
and  CLP  detector  performance. 
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Figure'*.  F t'l'oci  of \ai>  mg  N j- on  CLP  and  .VSl 
detector  pot formairce . 


V.  CONCLUSIONS 


I'oiir  ilcloctors  ofy  siiuisoiii  of  known  frcqucnfy  in  Ciaussian  while  noise  have  been 
analy/eil.  Two  of  these  ntili/.ed  tlie  Irequeney  response  of  an  adaptive  linear  predictive 
filler,  ami  two  employed  I'ourier  Iransfoiins  of  the  input  data.  The  results  ,diow  that  Ine 
CLP  proeessor  (I'ig.  3),  which  utilizes  the  deterministic  phase  of  the  l.PI'  frequency 
response,  offers  about  a l-dH  advantage  over  the  ASH  detector,  which  utilizes  averaged 
spectral  estimates  of  the  input  data  (Hig.  2).  The  ILP  detector  (I'ig.  4).  which  does  not 
utilize  the  LPI'  phase,  performs  about  the  same  as  the  ASH  detector.  The  ASH.  ILP.  and 
CLP  detectors  all  require  significantly  higher  SNK  than  the  optimal  detector,  which  employs 
a coherent  DP  T of  the  input  data.  In  all  cases,  the  CLP  and  ILP  adaptive  feedback  constants 
were  chosen  to  yield  the  best  possible  adaptive  filter  performance. 

The  work  presented  here  has  treated  the  use  of  adaptive  LPH  frequency  response  as 
a method  of  signal  detection.  Spectral  analysis  of  the  LPH'  output,  discussed  by  Zeidler  and 
Chabries  (Ref.  ‘)),  is  an  alternative  technique.  It  should  be  noted  that  the  detection  problem 
treated  in  this  paper  assumed  that  the  receiver  had  a considerable  amount  of  information 
concerning  both  the  signal  and  noise  characteristics.  As  suggested  by  Griffiths  (Ref.  10). 
conditions  of  greater  signal  uncertainty  (e.g..  nonstationary  noise  or  non-white  noise  of 
unknown  spectral  density)  may  result  in  adaptive  processor  perfonnance  superior  to  that 
of  nonadaptive  techniques. 
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APPl  NDIX:  LPI  WtStJU  l VI  ( TOK  MODI  1 

l itiiiiv  10  illiiNli.ilcs  llio  Imcai  pivilu  tui  iv>nMiU-ivil  id  this  ropoit.  A Ilioriitifili 
aisciiNsion  ol  llus  system  (also  kiuiwn  as  Ilio  ailaplive  line  enhaneei  1 has  been  piesenteil  by 
VVulrow  i ! ,jl  ( Kel.  1 1 aiul  l ieiehlei  ( Rcl.  1 1 ).  Only  a bne:  summary  is  piven  here.  1 he 
adaptive  lilter  ol  big.  10  is  a time  domain  digital  lilter  whieh  adr.ists  its  impulse  response  so 
as  to  mmnm/c  the  mean  square  ol  r , the  dilTerenee  between  the  desired  lilter  response, 

•nul  its  actual  vesponse.  y When  the  filter  ailinstment  is  complete,  y approximates  the 
minimum  mean  square  eris'r  prediction  of  x over  the  inter\al  ol  samples,  l-or  the  analysis 
which  follows,  it  is  assumed  that  x is  statistically  stationaiy  and  composed  of  a sine  wave 
m white  noise.  I'he  updating  algoritiun  for  the  k’**  iteration  ol  the  i’  ’ lilter  weight  is 


W(i,k  + n = \V(i,k)+  :MC(kl  xti,k) 


(1.0 


where  x(i,k)  i.v  Ihe  Jala  .sample  at  Ihe  i**’  filler  lap  during  the  k**'  iteration  and  W(i,k) 
denotes  the  i**'  weight  at  Ihe  k'*'  iteration.  Hie  constant  p controls  the  feedback  magnitude 
of  the  adaptive  system,  l or  stability,  p must  be  chosen  greater  than  zero  and  smaller  than 
the  reciprocal  of  Ihe  largest  eigenvalue  of  the  input  data  correlation  matrix,  R.  I he 
elements  of  R are 


inn 


= I'j.  1 x(m,k)  xtn,k)| , 0 < m,n  < N-1  , 


(14) 


where  R,„„  is  the  element  m Ihe  in’"  row  and  n’”  column,  and  1 1 • I denotes  expectation 

with  respect  to  Ihe  index  k.  Ihe  optimal  I IM  weight  vector,  resulting  in  the  minimum  mean 
square  oft,  is  given  by  Ihe  iliscrele  Wiener-llopf  equation; 


Wti) 


R 


-I 


ik  ’’k- 


()<i,k  < N-l  , 


(I.O 


where  IV  is  Ihe  k’*’  ciunponenl  of  1’,  Ihe  correl.ilion  vector  I'f  the  desired  filler  response  and 


Ihe  data  contained  in  the  filter,  i c , 


I’k  = I ,„| x(-fi,m I x(k,m )| 


N- 


I lo) 
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I'or  llio  case  of  a siiio  \va\o  of  luagiiitudo  A anil  frequoncN’  in  wliitc  noise  ot  total  powt*r 
o'  in  the  input  handwiilth.  the  elements  of  R are 


*^mn  ~‘’~^inn^ 


A- 


eos 


2;rfj^  tm 


-n)  A| 


()<  m.n  < N-1 , (17) 


wliere  A is  tlie  sampling  period  of  the  digital  LPF  system  and  6„,„  is  the  Kroneeker  della. 
The  eomponents  of  P are 


P 


m 


■) 

A- 


eos|  ^ttI,.  (m+6)  A1 


0<m<N-l.  (18) 


In  the  following  it  will  be  eonvenienl  to  assume  that  I.,  is  neither  /.ero  nor  the  Nyquist 
frequency  of  the  digital  system. 

Since  f^  is  known,  A and  the  predictor  weight  vector  length.  N,  can  be  chosen  to 
equal  an  integer  number  of  sine  wave  periods. 

fjjAN  = integer 

When  N is  chosen  in  this  way  R becomes  circulant  (Ref.  8).  Tlie  eigenvalues,  and 
eigenvectors,  lI^'’'\of  are  given  by 


j^(m)  = 


R()n  expl-j-»n«n/N| 


0<m  <N-1 


l,(m) 


I ^.-j2jrm/N  ^.-j27rm(N-l  )/N  | ^ 

l)<ni<N-l. 


(20) 


(21) 


where  T denotes  vector  transpose.  Note  that  and  ll''^"'”^ire  comple.x  conjugates. 

The  eigenvectors  are  complex  exponentials,  each  with  an  integral  number  of  cycles  in  the 
length  of  the  1 PI'  weight  vector.  Because  N has  been  appropriately  chosen  and  fj.  is  suitably 
restricted,  two  of  the  eigenvectors  correspond  to  frequency  f^.  1-quation  ( 1 5)  can  be  solved 
using  I qs.  ( 20)  and  (21).  I'he  optimal  weight  vector  is  found  to  be  a sinusoid  of  frequency  l\. 


W(i) 


SNR 


N 

1 + — SNR 


cosl 2jrfj^(i+6)  A) 


0<i<N-l  (22) 


where  SNR  denotes  the  sinusoid  to  noise  power  ratio.  Ihis  is  the  form  to  which  the  I PI 
weight  vector  would  converge,  given  an  infinitely  long  adaptation  interval  and  the  absence 
of  misadjustment  noise.  The  optimum  weight  vector  is  orthogonal  to  all  but  the  two  coniu- 
gate  eigenvectors  of  which  correspond  to  frequency  fj^  . 


The  eigenvalues  assoeiateci  with  these  eigenvectors  are  equal  and  are  given  by 


Treichler  ( Ref.  1 1 ) has  shown  tliat  these  eigenvalues  control  the  speed  with  which  tlie  LPI' 
weight  vector  converges  to  its  optimal  value.  They  are  inversely  proportional  to  the  con- 
vergence time  constant.  Assuming  the  weight  vector  to  be  initially  zero,  and  neglecting 
transient  behavior  as  the  predictor  initially  fills  with  data,  the  mean  value  of  the  i**’  weight 
at  the  k**'  update  is 

W(i.k)  = 11  - tl  -2^iX)^\  W(i).  (24) 

In  addition  to  this  mean  value,  the  weight  vector  also  contains  misadjustment  noise.  Widrow 
tRet's.  1,12)  has  shown  that  the  mean  value  of  this  noise  is  zero,  and  its  variance  is 
approximated  by  /ao*.  Misadjustment  is  uncorrelated  t'roin  weight  to  weight.  In  order  to 
extend  Widrow’s  analysis  of  detector  performance,  the  following  two  assumptions  will  be 
made  concerning  misadjustment  noise: 

( ( 1 ) Hie  mi.sadjustment  distribution  is  approximately  normal. 

(2)  Variance  of  the  misadjustment  is  equal  to  the  steady-state  value  of  at  all 
times,  and  is  not  intluenced  by  the  presence  or  absence  of  the  sinusoid. 

Reference  13  presents  data  from  e.xtensive  Monte  Carlo  e.xperiments  with  the  system 
of  Fig.  10.  These  show  that  the  steady  state  weight  misadjustment  is  normally  distributed 
when  the  filter  is  driven  by  uncorrelated  Gaussian  input  samples.  Tlie  results  presented  in  this 
paper  are  for  uncorrelated  Gaussian  input  plus  a sinusoid  of  less  than  -20  dB  relative  noise 
power.  The  presence  of  the  sinusoid  would  therefore  seem  to  be  unimportant  insofar  as 
misadjustment  is  concerned.  In  addition,  the  optimal  value  of  p (Eq.  10)  can  be  shown  to 
imply  that  the  observation  intend  of  N p input  data  samples  is  equal  to  about  1.2b  con- 
vergence time  constants  when  SNR  is  near  the  true  SNR.  Convergence  is  therefore  about 
70  percent  complete.  In  view  of  this,  the  use  of  the  steady  state  misadjustment  distribution 
ami  variance  a.ssumcd  here  would  seem  to  provide  a reasonable  representation  of  filter 
misadjustment. 

Under  the  above  assumptions,  the  noise  component  of  the  i**'  l.P!'  weight  is 

Pn  = N(’l0,  po“).  (2.*') 

'm 

I he  LI’F  weight  vector  is  modeled  by  the  sum  of  Fqs.  (24)  and  (25): 

W(i.k)  = 1 1 - (I  - 2pX)^|  W(i)  + Tjj  (2b) 

In  working  with  Fq.  (2b)  it  is  convenient  to  note  that,  for  the  situations  of  interest  in  this 
paper,  the  square-bracketed  term  is  of  the  form  | I - ( l-x)^  | . where  \ « 1 . y » 1 . and 
\*  y « I . This  term  is  closely  approximated  by  1 1 - e"’^^  I . as  can  be  shown  by  taking  the 
natural  logarithm  of  ( 1 - x)^  and  expanding  in  powers  of  x. 


SubNlitiiting  I q.  (22»  lor  \V(i)  aiul  evaluating  the  Dl  T of  f q.  (2M  at  fivqueiKy  reMili>  in 
j complex  number  with  ileterministic  aiul  random  parts,  (he  deterministic  part  is  due  to  the 
first  term  of  (-q.  (2M  and  is  equal  to 

N 

•7  SNR 

— ^ 1 1 - tl  exp  |j  :rr  f 6A|  . C’’) 

N 

1 + —SNR 

The  random  part  is  due  to  the  noise  component  of  I q.  (2o).  The  real  and  imaginary  compo- 
nents of  the  random  part  are  independent,  identically  distributed  normal  random  variables 

N , 

with  zero  mean  and  variance  equal  to  7 M**  ■ 


1.^ 
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